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Abstract 

We review some astrophysical and cosmological properties and implications of neutrino 
masses and mixing angles. These include: constraints based on the relic density of neu- 
trinos, limits on their masses and lifetimes, BBN limits on mass parameters, neutrinos 
and supernovae, and neutrinos and high energy cosmic rays. 
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1 Introduction 



The role of neutrinos in cosmology and astrophysics can not be understated [![]. They 
play a critical role in the physics of the early Universe, at temperatures scales of order 1 
MeV, and strongly determine the abundances of the light elements produced in big bang 
nucleosynthesis. They almost certainly play a key role in supernova explosions, and if they 
have mass, could easily contribute to the overall mass density of the Universe. At the present 
time, the only indicators of neutrino masses are from astrophysical sources, the inferred 
neutrino oscillations of neutrinos produced in the Sun, and those produced in cosmic-ray 
collisions in the atmosphere. Indeed, their elusive character has proven that a great deal of 
information on neutrino properties can be gained by studying their behavior in astrophsyical 
and cosmological environments. Here, we will try to elucidate some of these constraints on 
neutrino masses. 

In our discussion below, we will assume that the early Universe is well described by a 
standard Friedmann-Lemaitre-Robertson- Walker metric 

i]yl + r 2 (d6 2 + sin 2 6 dtf 



ds 1 = dt 2 - R 2 {t) 



1 — kr 2 
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We further assume that thermal equilibrium was established at some early epoch and that 
we can describe the radiation by a black body equation of state, p = p/3 at a temperature 
T. Solutions to Einstein's equations allow one to determine the expansion rate of the Uni- 
verse defined to be the Hubble parameter in terms of the energy density in radiation, the 
curvature and the cosmological constant. In the early Universe the latter two quantities can 
be neglected and we write 

^(f) 2 = 8?r 3* P (2) 

where the energy density is 

The present neutrino contribution to the total energy density, relative to the critical density 
(for a spatially flat Universe) is 

n* = - (4) 

Pc 

where p c = 1.06x 10~ 5 /i 2 GeV/cm 3 and h = i//100km/Mpc/s is the scaled Hubble parameter. 
For a recent review of standard big bang cosmology, see @. 



2 The Cosmological Relic Density of Stable Neutrinos 

The simplicity of the standard big bang model allows one to compute in a straightforward 
manner the relic density of any stable particle if that particle was once in thermal equilibrium 
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with the thermal radiation bath. At early times, neutrinos were kept in thermal equilibrium 
by their weak interactions with electrons and positrons. Equilibrium is achieved whenever 
some rate T is larger than the expansion rate of the Universe, or Tj > H . Recalling that 
the age of the Universe is determined by H' 1 , this condition is equivalent to requiring that 
on average, at least one interaction has occurred over the life-time of the Universe. On 
dimensional grounds, one can estimate the thermally averaged low-energy weak interaction 
scattering cross section 

M ~ g A T 2 /m* w (5) 

for T <C m\y Recalling that the number density scales as n oc we can compare the weak 
interaction rate T ~ n(av), with the expansion rate given by eqs. (fj) and (Q). Neutrinos 
will be in equilibrium when T w k > H or 

T 3 > yj8ir 3 N/90 m A w /M P (6) 

where M P = G~ N 1/2 = 1.22 x 10 19 GeV is the Planck mass. For iV = 43/4 (accounting for pho- 
tons, electrons, positrons and three neutrino flavors) we see that equilibrium is maintained 
at temperatures greater than 0(1) MeV (for a more accurate calculation see |J). 

The decoupling scale of 0(1) MeV has an important consequence on the final relic density 
of massive neutrinos. Neutrinos more massive than 1 MeV will begin to annihilate prior 
to decoupling, and while in equilibrium, their number density will become exponentially 
suppressed. Lighter neutrinos decouple as radiation on the other hand, and hence do not 
experience the suppression due to annihilation. Therefore, the calculations of the number 
density of light [m v < 1 MeV) and heavy [m v > 1 MeV) neutrinos differ substantially. 

The number of density of light neutrinos with m v < 1 MeV can be expressed at late 
times as 

p v = myVyii^ (7) 

where Y v = n u /n 1 is the density of i^'s relative to the density of photons, which today is 411 
photons per cm 3 . It is easy to show that in an adiabatically expanding universe Y v = 3/11. 
This suppression is a result of the e + e~ annihilation which occurs after neutrino decoupling 
and heats the photon bath relative to the neutrinos. In order to obtain an age of the Universe, 
t > 12 Gyr, one requires that the matter component is constrained by 

Qh 2 < 0.3. (8) 

From this one finds the strong constraint (upper bound) on Majorana neutrino masses: [[| 

m tot = ^m, < 28eV. (9) 

where the sum runs over neutrino mass eigenstates. The limit for Dirac neutrinos depends on 
the interactions of the right-handed states (see discussion below). As one can see, even very 
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Figure 1: Summary plot of the relic density of Dirac neutrinos (solid) including a possible 
neutrino asymmetry of r] u = 5 x 10~ n (dotted). 



small neutrino masses of order 1 eV, may contribute substantially to the overall relic density 
The limit (||) and the corresponding initial rise in Q u h 2 as a function of m v is displayed in 
the Figure [I] (the low mass end with m u < 1 MeV) . 

The calculation of the relic density for neutrinos more massive than ~ 1 MeV, is sub- 
stantially more involved. The relic density is now determined by the freeze-out of neutrino 
annihilations which occur at T < m u , after annihilations have begun to seriously reduce 
their number density ||. The annihilation rate is given by 

T ann = (av) ann n u ~ ^m u Tf' 2 e- m ^ T (10) 

where we have assumed, for example, that the annihilation cross section is dominated by 
uu — > ff via Z-boson exchange]] and (o~v) ann ~ m 2 /m A z . When the annihilation rate becomes 
slower than the expansion rate of the Universe the annihilations freeze out and the relative 
abundance of neutrinos becomes fixed. Roughly, Y v ~ {m(av) ann )~ 1 and hence Vt v h 2 ~ 
(o'v)ann~ 1 , so that parametrically Q u h 2 ~ l/m 2 ,. As a result, the constraint (§) now leads 
to a lower bound [^, ||, on the neutrino mass, of about m u > 3 — 7 GeV, depending on 
whether it is a Dirac or Majorana neutrino. This bound and the corresponding downward 
trend Q u h 2 ~ 1/m^ can again be seen in Figure III. The result of a more detailed calculation 
is shown in Figure |2] for the case of a Dirac neutrino. The two curves show the slight 
sensitivity on the temperature scale associated with the quark-hadron transition. The result 
for a Majorana mass neutrino is qualitatively similar. Indeed, any particle with roughly 
weak scale cross-sections will tend to give an interesting value of Qh 2 ~ 1. 

1 While this is approximately true for Dirac neutrinos, the annihilation cross section of Majorana neutrinos 
is p-w&ve suppressed and is proportional of the final state fermion masses rather than m v . 
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Figure 2: The relic density of heavy Dirac neutrinos due to annihilations 0. The curves are 
labeled by the assumed quark-hadron phase transition temperature in MeV. 



The deep drop in Q u h 2 , visible in Figure |l| at around m v = Mz/2, is due to a very strong 
annihilation cross section at Z-boson pole. For yet higher neutrino masses the Z-annihilation 
channel cross section drops as ~ \jm 2 vl leading to a brief period of an increasing trend in 
fl u h 2 . However, for m u > m w the cross section regains its parametric form (&v) ann ~ m 2 v 
due to the opening up of a new annihilation channel to H^-boson pairs ||, and the density 
drops again as Vl v h 2 ~ l/ra^. The tree level H^-channel cross section breaks the unitarity at 
around C(few) TeV |J however, and the full cross section must be bound by the unitarity 
limit (TtJ. This behaves again as 1/m 2 , whereby Q u h 2 has to start increasing again, until it 
becomes too large again at 200-400 TeV [II], |[| (or perhpas somewhat earlier as the weak 
interactions become strong at the unitarity breaking scale). 



3 Neutrinos as Dark Matter 

Based on the leptonic and invisible width of the Z boson, experiments at LEP have de- 



termined that the number of neutrinos is N u = 2.9841 ± 0.0083 |p| . Conversely, any new 
physics must fit within these brackets, and thus LEP excludes additional neutrinos (with 
standard weak interactions) with masses m v < 45 GeV. Combined with the limits displayed 
in Figures [3] and we see that the mass density of ordinary heavy neutrinos is bound to be 
very small, Q u h 2 < 0.001 for masses m u > 45 GeV up to m u ~ (9(100) TeV. 

A bound on neutrino masses even stonger than Eqn. (^) can be obtained from the recent 
observations of active-active mixing in both solar- and atmospheric neutrino experiments. 
The inferred evidence for — v T and v e — T mixings are on the scales m 2 ~ 1 — 10 x 10~ 5 
and ml ~ 2-5 x 10~ 3 . When combined with the upper bound on the electon-like neutrino 



mass m u < 2.8 eV [12|, and the LEP-limit on the number of neutrino species, one finds the 
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constraint on the sum of neutrino masses: 



0.05 eV < m tot < 8.4 eV. (11) 

Conversely, the experimental and observational data then implies that the cosmological 
energy density of all light, weakly interacting neutrinos can be restricted to the range 

0.0005 < tt u h 2 < 0.09. (12) 

Interestingly there is now also a lower bound due to the fact that at least one of the neutrino 
masses has to be larger than the scale m 2 ~ 10~ 3 eV 2 set by the atmospheric neutrino 
data. Combined with the results on relic mass density of neutrinos and the LEP limits, the 
bound ( |12"D implies that the ordinary weakly interacting neutrinos, once the standard dark 
matter candidate [13], can be ruled out completely as a dominant component of the dark 



matter. 

This conclusion can be evaded if neutrinos are Dirac particles, and have a nonzero asym- 
metry however, since then the relic density could be governed by the asymmetry rather than 
by the annihilation cross section. Indeed, it is easy to see that the neutrino mass density 
corresponding to the asymmetry r] u = {n u — n^/n^ is given by |L4 



p = m v T] u n^, (13) 

which implies 

Vt u h 2 ~ 0.004 r) vW (m^/GeV) . (14) 

where t] v \q = 10 10 ^. We have shown the behaviour of the energy density of neutrinos with 
an asymmetry by the dotted line in the Figure |l|. At low m„, the mass density is dominated 
by the symmetric, relic abundance of both neutrinos and antineutrinos which have already 
frozen out. At higher values of m v , the annihilations suppress the symmetric part of the 
relic density until Q u h 2 eventually becomes dominated by the linearly increasing asymmetric 
contribution. In the figure, we have assumed an asymmetry of r] u ~ 5 x 10~ n for neutrinos 
with standard weak interaction strength. In this case, Q u h 2 begins to rise when m v > 20 
GeV. Obviously, the bound (§) is saturated for m„ = 75GeV/r] uW . 

There are also other cosmolgical settings that give rise to interesting mass constraints 
on the eV scale. Indeed, light neutrinos were problematic in cosmology long before the 
imporoved mass limits leading to ( |T2"D were established, due to their effect on structure 
formation. Light particles which are still relativistic at the time of matter domination erase 
primordial perturbations due to free streaming out to very large scales []T3]| . Given a neutrino 



with mass m„, the smallest surviving non-linear structures are determined by the Jean's mass 

Mj = 3 x io 18 4rk • ( 15 ) 
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Thus, for eV mass neutrinos the large scale structures, including filaments and voids ]TJj, |T7 
must form first and galaxies whose typical mass scale is ~ 1O 12 M are expected to fragment 
out later. Particles with this property are termed hot dark matter (HDM). It seemed that 
neutrinos were ruled out because they tend to produce too much large scale structure 



and galaxies formed too late flT], jl9| , at z < 1, whereas quasars and galaxies are seen out to 
redshifts z > 6. 

Subsequent to the demise of the HDM scenario, there was a brief revival for neutrino dark 
matter as part of a mixed dark matter model, using now more conventional cold dark matter 
along with a small component of hot (neutrino) dark matter. The motivation for doing this 



was to recover some of the lost power on large scales that is absent in CDM models [20 



However, galaxies still form late in these models, and more importantly, almost all evidence 
now points away from models with Q m = 1, and strongly favor models with a cosmological 
constant (ACDM). 

Combining the rapidly improving data on key cosmological parameters with the better 
statistics from large redshift surveys has made it possible to go a step forward along this 
path. It is now possible to set stringent limits on the light neutrino mass density Q u h 2 , and 
hence on neutrino mass based on the power spectrum of the Ly a forest ||21|| , m tot < 5.5 
eV, and the limit is even stronger if the total matter density, Q m is less than 0.5. This 
limit has recently been improved by the 2dF Galaxy redshift |22| survey by comparing the 
derived power spectrum of fluctuations with structure formation models. Focussing on the 
the presently favoured ACDM model, the neutrino mass bound becomes m tot < 1.8 eV for 
Vt m < 0.5. 

Finally, right handed or sterile neutrinos may also contribute to the dark matter. The 
mass limits for neutrinos with less than full weak strength interactions are relaxed [23]. For 
Dirac neutrinos, the upper limit varies between 100 - 200 eV depending on the strength 
of their interactions. For Majorana neutrinos, the limit is further relaxed to 200 - 2000 
eV. This relaxation is primarily due to the dilution of the number density of super-weakly 
interacting neutrinos due to entropy production by decay and annihilation of massive states 
after their decoupling from equilibrium ||24j| . Such neutrinos make excellent warm dark 
matter candidates, albeit the viable mass range for galaxy formation is quite restricted PS] . 



4 Neutrinos and Big Bang Nucleosynthesis 

Big bang nucleosynthesis is the cosmological theory of the origin of the light element isotopes 
D, 3 He, 4 He, and 7 Li ||26|| . The success of the theory when compared to the observational 
determinations of the light elements allows one to place strong constraints on the physics 
of the early Universe at a time scale of 1-100 seconds after the big bang. 4 He is the most 
sensitive probe of deviations from the standard model and its abundance is determined 
primarily by the neutron to proton ratio when nucleosynthesis begins at a temperature of 
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~ 100 keV (to a good approximation all neutrons are then bound to form 4 He). The ratio 
n/p is determined by the competition between the weak interaction rates which interconvert 
neutrons and protons, 



p + e <->• n + v e , n + e + <-> p + v e , n^p + e +z/ e (16) 
and the expansion rate, and is largely given by the Boltzmann factor 

n/p ~ e - {m "- mp)/ ^ (17) 

where m n — m p is the neutron to proton mass difference. As in the case of neutinos discussed 
above, these weak interactions also freeze out at a temperature of roughly 1 MeV when 



Gf Tf ~ r wk (7» = H(T f ) ~ JG N NT/ (It 



The freeze-out condition implies the scaling Tj ^ yN . From Eqs. ( |i7j ) and fllq) , it is then 
clear that changes in N, caused for example by a change in the number of light neutrinos N u , 
would directly influence n/p, and hence the 4 He abdundance. The dependence of the light 
element abundances on N u is shown in Figure |3| [27| , where plotted is the mass fraction of 
4 He, Y, and the abundances by number of the D, 3 He, and 7 Li as a function of the baryon- 
to-photon ratio, rj, for values of N u = 2 — 7. As one can see, an upper limit to Y, combined 



with a lower limit to n will yield an upper limit to N v f28[ . 

Assuming no new physics at low energies, the value of n is the sole input parameter 
to BBN calculations. It is fixed by the comparison between BBN predictions and the ob- 



servational determinations of the isotopic abundances [p9j] . From 4 He and 7 Li, one finds 



a relatively low value [^(J, [29j for n ~ 2.4 x 10~ 10 corresponding to a low baryon density 
VLsh 2 = 0.009 with a 95% CL range of 0.006 - 0.017. Deuterium, on the other hand, implies 
a large value of rj and hence a large baryon density: n ~ 5.8 x 10~ 10 and VLbIi 2 ~ 0.021 
with a 95% CL range of 0.018 - 0.027. The value of the baryon density has also been 
determined recently from measurements of microwave background anisotropics. The recent 
result from DASI |fi] indicates that Q B h 2 = 0.022±g;gg|, while that of BOOMERanG-98 [g2 



Q B h 2 = 0.02lljj:oo3 (using la errors). 

With the value of n fixed, one can use He abundance measurements to set limits on 
new physics. In particular one can set upper limits the number of neutrino flavors. Taking 
Y p = 0.238 ± 0.002 ± 0.005 (see e.g. [0]), we show in Figure ^ the likelihood functions for 



N v based on both the low and high values of n pTl. The curves show the impact of an 



increasingly accurate determination of r\ from 30% to 3%. If one assumes a 20% uncertainty 
in rj (the current uncertainty level), these calculations provide upper limits of 

N v < 3.9 rj = 2.4 x 10~ 10 

N v < 3.6 n = 5.8 x lO" 10 (19) 
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Figure 3: The light element abundances as a function of the baryon-to-photon ratio for 
different values of N v p7 |. 



at the 95% CL. Although, as noted above, LEP has already placed very stringent limit to 
N u , the limit flT9]) is useful, because it actually applies to the total number of new particle 
degrees of freedom and is not tied specifically to neutrinos. In fact, more generally, the 
neutrino limit can be translated into a limit on the expansion rate of the Universe at the 
time of BBN, which can be applied to a host of other constraints on particle properties. 



4.1 BBN limits on Neutrino masses and lifetimes 

As discussed above, the nucleosynthesis prediction for light element abundances is sensitive 
to the changes in the expansion rate of the universe, which depends on the energy density of 
the universe during the BBN era (0). This extra energy density could be in the form of new 
massless degrees of freedom, in which case their number is directly constrained by equation 
(|19|) . Equally well the extra energy density could reside in the form of massive long lived but 
unstable neutrinos, in which case nucleosynthesis provides interesting constraints on their 
masses and life-times. 

We already pointed out that the relic density of neutrinos strongly depends on whether 
they decouple while relativistic or nonrelativistic. Here the calculations also depend on 
how the neutrino life-times relate to the BBN time-scale of about 100 seconds. Also, in 
order to get reliable results for the light element abundances, one must keep track of the 
induced perturbations (electron neutrino heating) in the weak reaction rates ( jl6|) in addition 
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Figure 4: (a) The distribution in N v assuming a value of rj 
and the CBI measurement of the microwave background anisotropy [27 



2.4 x 1(T 10 from 4 He and 7 Li 
The curves show 



the effect of the expected increased accuracy in the CMB determination of rj. (b) (b) As 
in (a), but assuming a value of rj = 5.8 x 10~ 10 from D and the DASI and BOOMERanG 
measurements of the microwave background anisotropy. 



to computing changes in the expansion rate. Nevertheless, even in this case it is customary 
to measure the change in helium abundance in units of equivalent effective neutrino degrees 
of freedom N eS , such that the limit (fity) can be applied on N efi (AY(m u , r u )). 

When the neutrino life-time is much larger than 100 seconds, neutrinos are effectively 
stable on a nucleosynthesis scale |34|, [3^, While accounting for changes in the rates in 



Eq. (]I6|) is important for the detailed bounds, the bulk behaviour of N eS (m u ) is dictated 
by the neutrino mass contribution to the energy density. Obviously, when m u < 0.1 MeV, 
neutrinos are effectively massless during BBN, and N e s(m u ) — > 3 when m v — > 0. For masses 
in excess of 0.1 MeV, but below the neutrino decoupling temperature of 0(iew) MeV, their 
number density is unsuppressed and their mass density can be large during BBN, causing 
N c q to increase. For yet larger masses however, the Boltzmann factor shown in (|T0l ) begins 
to suppress the mass density and eventually turns N e g down again. This behaviour is shown 
in Figure [5] for a massive Dirac and Majorana type tau-neutrino ||35|| . The bound (|19D yields 
an excluded region for stable neutrino masses centered around (9(few) MeV. For N v < 3.6 
the lower bound is m u > 42 MeV (Majorana) and m u > 30 MeV (Dirac) ||35|| . This is only 
relevant for u T , and is complementary to the present laboratory limit on the r-like neutrino, 



m u < 18 MeV [38|. Due to contributions from pion decays and inverse decays to neutrinos, 
the upper bound from BBN depends on the QCD-phase transition temperature, Tqcd, and 
is also different for r- and /^-neutrinos because of their different scattering rates off muons. 



Imposing again the constraint N v < 3.6, and taking Tqcd = 200 MeV gives p9fl : 

m v < 230 KeV // - like 
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Figure 5: Plot of the effective number of neutrino degrees of freedom during BBN for a Dirac 
(dashed) and a Majorana neutrino (solid) [ |35|] . 



m u < 290 KeV r - like. (20) 

The laboratory limit on the muon-like neutrino, for comparision, is m u < 170 keV j|(J. To 
improve this, the BBN limit should be improved to N„ < 3.4 |39| . 



When the neutrino life-time is small or comparable to the nucleosynthesis time scale, one 
has to account for neutrino decay processes as well. This involves solving for the distributions 
of the final state decay products, which might include new particles like majorons, and their 
possible direct effect on BBN. (For example energetic photons would cause the dissociation of 
the newly generated light nuclei.) Such calculations have been done by many groups jy], 42] , 



and the results are given in exclusion plots in the mass-vs-life-time plane. Constraints are 
possible for masses of order 0(1) MeV and life-times of order 0(1) second. In Figure we 
show a constraint on tau-neutrino masses and life-times as an example (data taken from the 



reference |42|). 



4.2 BBN limits on Neutrino mixing parameters 

Despite losing the competitive edge w.r.t. masses and life-times, BBN continues to put inter- 
esting limits on other neutrino mass parameters, relevant for neutrino oscillations. Indeed, 
the LEP-limit of course only applies to neutrinos with weak interactions, while neutrinos 
without weak interactions, or sterile neutrinos, have been proposed in many different con- 
texts over the years. At present, a prime motivation to introduce sterile neutrinos is to 
explain the LSND neutrino anomaly |Q in conjunction with the solar and atmospheric 



neutrino deficits. 

BBN on the other hand is sensitive to any type of energy density changing the expan- 
sion rate in the (9(0.1-1) MeV range, irrespective of their interactions. It is therefore very 
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Figure 6: Plot of BBN-constraint on masses and life-times of an unstable tau neutrino. 
Contours are labeled by one (dashed) and two sigma (solid) deviations from the observed 
value of Y p . The upper right corner is excluded due to too much and the lower part of the 
graph by too little 4 He being produced. 

interesting to observe that even if no sterile neutrinos were created at very early times, they 
could be excited by mixing effects in the early universe prior to nucleosynthesis. The basic 
mechanism is very simple. Suppose that an active state u a (a = e, //, r) mixes with a sterile 
state v s . That is, the neutrino mass matrix and hence the Hamiltonian, is not diagonal in 
the interaction bases. The mixing is further affected by the forward scattering interactions 
with the background plasma, felt by the active state. As a result, even though a neutrino 
state was initially produced in purely active projection, after some time t it has become some 
coherent linear combination of both active and sterile states: 



The coherent evolution of this state is interrupted by collisions, which effect a sequence of 
quantum mechanical measurements of the flavour content of the propagating state. Since 
the sterile state has no interactions, each measurement is complete, and collapses the wave- 
function to the sterile state with a probability P Ve ^ Us {t) = |c s (t)| 2 . As a result, the sterile 
states are populated roughly with an average rate 



where T Va is the weak interaction rate of the active state z/ a , and we have assumed that 
the oscillation time is short in comparision with the collision time scale. The matter mixing 



c e (t)u e + c B (t)u s . 



(21) 




(22) 
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Figure 7: Plotted are the BBN constraints on the active sterile neutrino mixing parame- 
ters jUJ for v a = v e (left) and u a = v^ T (right). Regions to the right from the contours 
labeled by the bound on effective neutrino degrees of freedom 5N eS = N cS — 3 are excluded. 
Shown are also the current regions corresponding to active-sterile mixing parameters for 
atmospheric (ATM) and the large mixing angle (LMA) solar neutrino solutions. 



angle 8 m is given by f|4 



9 „ sin 2 29n , . 

sin 2 26 m = — 23 

1 — 2x cos 20 + x 

where sin2#o is the mixing angle in vacuum and \ = 2p\V\/5m 2 where 5m 2 is the mass 
squared difference of the vacuum mass eigenstates, p ~ T is the momentum and \V\ is 
the matter induced effective potential to the Hamiltonian |5tj, The weak rate scales 
as T Ua ~ T 5 . Moreover % ~ T 6 at very high temperatures, which causes a strong matter 
suppression for mixing and hence r„ s ~ T~ 7 . At very small temperatures 6 m — > 6 on the 
other hand, and hence Y Vs ~ T 5 . The rate is thus suppressed both at very large and at very 
small temperatures \ In the intermediate region of O(few) MeV however, T VfS can exceed 
the expansion rate bringing a significant amount of sterile neutrinos into equilibrium. An 
accurate treatment of the problem requires a numerical solution of the appropriate quantum 
kinetic equations, and the results depend on whether the mostly active state is heavier 
(5m 2 < 0) or lighter (5m 2 > 0) of the mixing states. We show the results of such a 
calculation in Figure |7| below pifl . The lines are labeled by constant effective number of 
degrees of freedom during BBN: 5N V = N u — 3. The most recent limits corresponding to 
(P~5|) can be interpolated from the curves shown. It is the area above the curves which is 
excluded by BBN-limit. 

'in 



The BBN-limit can be converted to an upper bound on the sterile neutrino flux [47 
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the atmospheric and the solar neutrino observations. Using N v < 3.6 one finds 

sin 2 6^ s < 0.03 (Atmospheric) 

sin 2 es <0.06 (Solar LMA), (24) 

whereas the bounds from the atmospheric and solar neutrino experiments are about an order 
of mangitude weaker: sin 2 6^ s < 0.48 and sin 2 6 es < 0.72. 

The constraints shown in the Figure |7| and in Equation ([24]) depend on the assumption 
that the primordial lepton asymmetry is not anomalously large In ref. [[49] it was sug- 



gested that a large effective asymmetry violating this assumption could actually be generated 
by oscillations given a particular neutrino mass and mixing hierarchy. For a while these ideas 
generated a lot of interest, as they would have allowed reconciling all observed anomalies 
(including LSND) with the nucleosynthesis constraints. However, in these scenarios at least 
one of the active states would have to be much heavier than the two others, which is not 
allowed by the atmospheric and solar neutrino flux observations. As a result, the bounds 
(f24]) hold, and in particular nucleosynthesis is very much at odds with the possible existence 
of the LSND-type sterile state. To see this, observe that creating a large enough effective 
mixing between and z/ e to explain the anomaly |43| would require a sterile intermediare 
with m Us ~ 1 eV, and sin 2 26^ e ~ \ sin 2 2 6^ sin 2 26 sc > 10~ 2 (5m 2 /eV 2 )" 2 . In other words at 
least one of the active-sterile mixings should satisfy sin 2 26( P:e ) s > 0.15(5m 2 /eV 2 ) -1 , which 
is well within the BBN excluded regions shown in Figure [7]. (It would be excluded even by 
Acfr ^, 3.9, although we do not show that contour in Figure |7|). 

It should finally be noted that active-active type oscillations have hardly any effect on the 



expansion rate or the weak interaction rates |p0 |, and hence are not constrained in the above 



sense by BBN. However, large mixing angle active-active oscillations could equilibrate the 
lepton asymmetries prior to BBN. This has been used to put strong bounds on muon- and 



tau- lepton asymmetries [51], which exclude the possibility of the degenerate nucleosynthesis. 



5 Neutrinos and Supernovae 

Neutrinos have for long been known to play important role in the physics of supernovae. To 
be sure, it is clear that by far the largest part, roughly 99%, of the gravitational binding 
energy of about 3 x 10 53 erg involved in the explosion of a type II supernova is carried out 
by neutrinos, while just 1% powers the shock wave responsible for blowing out the mantle 
of the star, and only a tiny fraction of about 0.1% escapes in the form of light responsible 
for the spectacular sights observed in the telescopes watching the sky. 

The formation of the neutrino burst in the collapse of a type II supernova is rather well 
understood. The temporal structure of the burst and the energy spectrum of the emitted 



neutrinos can be computed fairly well [p2|| . Existing or planned large scale neutrino detec- 



tors |53[ can be used to observe deviations from these signatures and to obtain interesting 
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Figure 8: Time evolution of neutrino luminosities and average energies: v x represents the 
spectrum of u T , and v T . Figure taken from |53 . 



information on neutrino masses and mixing parameters |p4| , given a future observation of a 
Galactic supernova. We show an example of a compilation of neutrino fluxes and spectra in 
Figure ||. 

A number of constraints on new physics and in particular on neutrino parameters have 
already been deduced from the famous SN1987A event in the Small Magellanic Cloud. Of 
these, perhaps the most direct is the upper bound on the neutrino mass derivable from 
the maximum duration of the observed duration of the neutrino pulse of about 10 seconds. 
Given the initial energy spectrum of neutrinos and the distance to the supernova, one can 
compute the expected spread in the arrival times of the neutrinos to earth as a function of 
neutrino mass. Comparing the predictions with the observations has been shown to yield 
the bound [pSfl 

m Ue <Q- 20 eV. (25) 

The observed pulse length also lends to the classic cooling argument: any new physics that 
would enhance the neutrino diffusion such that the cooling time drops below the observed 
duration, must be excluded. Cooling arguments have been used to set limits on various 
neutrino properties |5B| such as active-sterile neutrino mixing |57| and neutrino magnetic 
moments. The magnetic dipole moment bound in particular was recently revised by Ayala 



etal. 58 to 



< 1 - 4 x 10 



-12, 



(26) 



which is a couple of orders of magnitude more stringent than the best laboratory bounds, and 
comparable to the bound coming from the globular cluster red giant cooling arguments pi . 
which give \i Ve < 3 x 10~ 12 /j jB . 

At present, most of the activity concerning neutrinos in supernovae has focussed on 
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the effects of neutrino transport in supernova explosion dynamics rather than with finding 
constraints on neutrino mixing parameters or masses. Indeed, the details of the physics 
responsible for the actual visibly observed supernova explosion, including blowing out the 
stellar mantle, are not very well understood. In particular, the shock wave, which forms deep 
within the iron core as the infall of matter is reversed due to the stiffening of the nuclear 
matter equation of state, is typically found to be too weak to explode the star. This is 
believed to be due to energy loss from the shock and dissociating iron nuclei on its way out 
from the core to the mantle. In the popular "delayed explosion scenario" , the stalling shock 
wave is rejuvenated by energy transfer to the shock from the huge energy-flux of neutrinos 
free streaming away from the core. Recent numerical simulations including diffusive neutrino 
transport do not verify this expectation however [60]; while neutrinos definitely help, they 
do not appear to solve the problem. These results are not conclusive because the diffusive 
transport equations used so far [BO] did not include all relevant neutrino interactions, most 



notably the nuclear brehmsstralung processes [pi] , p2| . Furthermore, processes other than 



diffusive processes, such as convective flows in the core and behind the shock, appear to play 
an important role as well [|6^| . 

It is of course possible that a succesfull SN explosion requires help of some new physics 
to channel energy more efficiently to the shock, and neutrino-oscillations have already been 
considered for this role ||63|| . The idea is that z/ M and z/ r , interact more weakly and hence 
escape more energetic from deeper in the core than do electron neutrinos. Arranging mixing 
parameters so that T turn resonantly to v e in the mantle behind the shock could increase 
the energy deposited to the shock significantly. Unfortunately the mass difference needed 
for the resonant transition would be very big: 

5m 2 ~ 2EV eS ~ 1.5 x 10 5 p e ^°° eV 2 (27) 

where p e is the electron density in units 10 10 g/cm 3 , which at the shock front is about 10~ 3 , 
and El°° is the neutrino energy in units of 100 MeV. So one gets m v > 10 eV, which is 
excluded by the present data. A similar idea was behind the suggestion [Bl|, that neutrino 
magnetic moments induce resonant transitions from vr (which escape energetic from deep 
in the core) to vl behind the mantle. This mechanism could actually be used to evade the 
bound (p6|) , but as it demands a magnetic moment of order ji v ~ lO -11 ^ it has problems 



in coping with the red giant cooling bound |59]. 

Sterile neutrinos could also be relevant for supernovae by alleviating the problems with 
the r-process nucleosynthesis, which is thought to be responsible for creation of the most 
heavy elements. In the standard SN calculations the r-process nucleosynthesis is not effective 
due to too efficient de-neutronization by the processes v e + n — > e~ +p. If electron neutrinos 
mixed with a sterile state however, these processes could be made less effective, increasing 
the neutron density in the mantle, and hence improving the r-process efficiency [ |35|| . 

Finally, there is also the old problem of the "kick" -velocities of pulsars (neutron star 
remnants of supernova explosions). It has proven difficult to arrange for these velocities, 
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which average around 450 km/sec, based on the normal fluid dynamics in asymmetrical 
SN-explosions. The momentum carried out by neutrinos p u ~ E u , on the other hand, is 
about 100 times larger than the pulsar kinetic energy, so that a mere one per cent asym- 
metry in the neutrino emission would be enough to power the pulsar velocities. Interesting 
attempts have been made to explain such asymmetric emission by an asymmetric distribu- 
tion of inhomogeneities in the SN magnetic fields, combined with a large neutrino magnetic 



moment [64], or just the magnetic field induced deformation of the neutrino spheres [56 



However, the former would again need probably too large magnetic moment to work, and a 
detailed analysis of the latter suggests that the a symmetric flux is very suppressed, requiring 
perhpas unrealistically large magnetic fields; according to [67] the field needs to be in excess 
of 10 17 G, while [^8| argue that a field of 10 14_15 G would suffice. The true nature of the 
physics explaining the kick velocities may remain ambiguous for some time to come, but a 
neutrino solution looks definitely appealing from the pure energetics point of view. 

Before we conclude this section, we would like to mention, several other astrophysical 
limits on neutrino properties. A sure limit to the mean life/mass ratio is obtained from 
solar x- and 7- ray fluxes [69]. The limit is r/m Vl > 7 x 10 9 s/eV. This is far superior 
to the laboratory bound of 300 s/eV |fr0| . Other much stronger limits (> O(10 15 ) s/eV) 
are available from the lack of observation of 7-rays in coincidence with neutrinos from SN 



1987A ||71|| . This latter limit applies to the heavier neutrino mass eigenstates, z/2 and V3, as 
well. 



6 Neutrinos and Cosmic Rays 

One of the most interesting puzzles in astrophysics today concerns the observations of ultra 
high energy cosmic rays (UHECR), beyond the so called Greisen-Zatsepin-Kuzmin (GZK) 
cutoff 

Eqzk — 5 x 10 19 eV. (28) 

The problem is that cosmic rays at these energies necessarily need to be of extragalactic 
origin, since their gyromagnetic radius within the galactic magnetic field far exceeds galactic 
dimensions. Yet, the attenuation lengths of both protons and photons are rather small in 
comparision with intergalactic distances, and neither can have originated by further than 
about 50 Mpc away from us, due to their scattering off the intergalactic cosmic photon 
background. As a result, one would expect that the cosmic ray spectrum would abruptly 
end around E ~ Eqzk due to scattering off of the microwave background. This cutoff is 
represented by the dotted line in Figure ^. In contrast, several groups, most notably the 
HiRes |7_2j and AGASA J73| collaborations, have reported events with energies well above 
the GZK-cutoff: the latest compilation of AGASA, for example, contains 10 events above 
the scale E > 10 20 eV, observed since 1993. 
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Figure 9: Seeled by E 3 spectrum of highest energy cosmic rays near the GZK cutoff. The 
dotted line corresponds to the expectation from uniformly distributed extragalactic sources, 
and the solid line shows the prediction of a Z-burst model of ref . ]74| . (Figure modified from 
the orginal found at the AGASA web page |73| .) 



The origin of UHECR's has been the subject of lively discussions over the last few years. 
Apart from having astrophysical orgins, being acclerated in extreme environments at extra- 
galactic distances (in AGN's, GRB's or Blazars), they have been attributed to the decay 
products of very heavy particles or of topological defects. All these explanations have prob- 
lems, however. Astrophysical explanations face the difficult task of accelerating particles to 
the extreme energies required, with little or no associated sub-TeV-scale photonic compo- 
nent (as none have been ever observed). This is in addition to the above mentioned problem 
of the propagation of UHECR's over extragalactic distances. Decay explanations are some- 
what disfavoured by the growing evidence of doublets and triplets in AGASA data, and by 
the correlation between the UHECR arrival directions with far compact Blazars [75], which 
appear rather to be pointing towards astrophysical origin. 

The attenuation problem for extragalactical UHECR's can be avoided however, if they 
cross the universe in the form of a neutrino beam, since neutrinos travel practically free over 
super-Hubble distances. Indeed, should this be the case, the initial ^uhe's could occasionally 
interact with the cosmological relic neutrino background close to us, giving rise to "Z-bursts" 
of hadrons and photons [76 1, which then would give rise to the observed UHECR events. 
Indeed, given a neutrino mass m u , such a collision has sufficient CM-energy for resonant Z 
production if 



Ev = M*. ~ 4.2 x 10 21 eV ( —) . (29) 
2m u \m u J 

The requirement of super-GZK-energies for the initial i^uhe beam leads immediately to the 
interesting mass scale for neutrinos: m v ~ 0(1) eV. 
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Z-burst models have been extensively studied lately [|77|| . For example, in ref. [74] it was 
shown that a Z-burst model with m v = 0.07 eV, corresponding to a degenerate neutrino 
spectrum, could reproduce the AGASA-data including the spectral features, such as the 
"ankle" and the "bump" observed at E < Egzk- The best fit model of ref. [74] is shown 
by the solid line in Figure |9[ While the agreement with AGASA data is good, it should 
be noted that this model predicts that cosmic ray primaries are exclusively photons above 
E > 10 20 eV, whereas the E ~ 3 x 10 20 eV event observed by Fly's Eye is almost certainly 
not caused by a photon . The model also predicts a large increase of the cosmic ray flux 
above fewxlO 20 eV (as a direct result of the huge initial energy needed: E UuiiE ~ 6 x 10 22 
eV) , which excacerbates the already outstanding problem of the origin of UHECR's. These 
problems could be ameliorated by assuming somewhat larger neutrino masses, and it has 
been argued by Fodor etal. ]77| , that the Z-burst scenario can already be used to constrain 
the mass, plausibly to within m v ~ 0.08 — 1.3 eV, in very good agreement with other mass 
determinations. The analysis of ref. |77| , was restricted to Z- resonance interactions however, 
while for higher CM-energies the cross section || for pair production of gauge bosons vv — > 
ZZ, WW becomes important. These reactions have been shown to give accptable solutions 
with larger neutrino masses m v > 3 eV [[79 1. 

In summary, although the origin of UHECR's remains a mystery today, it is almost certain 
that if it has to do with extragalactic sources, neutrino physics plays a very important role 
in the interpretation of these events. 
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